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Abstract. Crystalline acetanilide is a an organic solid with peptide bond structure similar 
to that of proteins. Two states appear in the amide I spectral region having drastically differ- 
ent properties: one is strongly temperature dependent and disappears at high temperatures 
while the other is stable at all temperatures. Experimental and theoretical work over the 
past twenty five years has assigned the former to a selftrapped state while the latter to an 
extended free exciton state. In this article we review the experimental and theoretical devel- 
opments on acetanilide paying particular attention to issues that are still pending. Although 
the interpretation of the states is experimentally sound, we find that specific theoretical 
comprehension is still lacking. Among the issues that that appear not well understood is the 
effective dimensionality of the selftrapped polaron and free exciton states. 

1 Introduction 

The first numerical experiment performed by Fermi, Pasta and Ulam in 1955 proved to be the beginning 
of an exciting, non-reductionist branch of modern science focusing on nonlinearity and complexity in a 
variety of physical systems 1 1 1. In this path, the discovery of the soli ton in the sixties by Zabusky and 
Kruskal and the understanding of many of its mathematical properties led in the seventies in the explo- 
ration for physical applications. Proposals for soliton modes in Josephson junctions, optical fibers and 
other physical systems where suggested at this time, while Davydov introduced the novel idea that soli- 
tons may have a direct impact to biology as well |3 4 5 1. More specifically Davydov proposed that the 
energy released in the process of ATP hydrolysis becomes vibrationally selftrapped and forms a soliton. 
The main feature that made solitons appealing to biology is their property of dispersiveless transport; if 
a soliton-like packet of energy forms in a macromolecule, it may propagate without losses and enable 
long range coherent energy transfer This feature could be critical in the internal energetics of proteins 
where energy deposited on given sites produces large effects at relatively distant locations. Although 
soliton propagation experiments in proteins were not possible at this time, Careri and coworkers were 
investigating independently crystalline acetanilide (ACN), a molecular solid that has peptide bonds and 
a structure similar to that of a protein |6 1 (Fig 1). They found in infrared absorption experiments that an 
"anomaly" appears in the amide-I spectral region manifested through a strongly temperature dependent 
spectral peak |7 8| (Fig. 2). In the early eighties Scott and collaborators surmised that this peak was 
related to a Davydov soliton M7I9I10I . Davydov's original soliton idea coupled to the experimental find- 
ings of Careri and theoretical picture of Scott for ACN lead to a funy of activity during the rest of the 
decade as well as in the early nineties. Brown et al. examined critically the general "Davydov sohton" 
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hypothesis 01 1I12II while in acetanilide Scott as well Alexander and Krumhansl soon realized that it 
was a small polaron rather than an extended soliton responsible for the anomalous peak II13I14I15I10I . 
The small polaron picture corroborated the neutron and absorption experiments of Barthes and col- 
laborators. II16I17I . During the same period, explicit lattice discreetness was introduced in nonlinear 
approaches and, as a result, a new entity was introduced, viz. that of an intrinsic localized mode or 
discrete breather 1 18 1. The latter modes where seen in many instances to describe better nonlinear local- 
ization than continuous solitons. The over thirty-year long history of the fertile Davydov soliton idea is 
landmarked presently by a new generation of experiments; the pump-probe experiments of Hamm and 
collaborators have not only corroborated the acetanilide picture but also produced the first experimental 
sign for nonlinearly localized states in proteins 1191201211221 . In the present work we will attempt to 
summarize the basic theoretical and experimental steps in the Davydov soliton idea as related specif- 
ically to acetanilide. We prefer to focus only on the latter since the excellent review by Scott in 1992 
covers much of the generalities on the theoretical Davydov soliton idea, while, additionally, acetanilide 
provides a tangible system where both theory and experiments have been employed. 

Nonlinear excitations in biomolecules are typically associated with with objects such as solitons, 
polarons or discrete breathers (DBs). All three modes are similar in nature in that they arise due to the 
actual or effective presence of some type of nonlinearity in the equations of motion. In the original 
formulation of the Davydov problem, a vibrational excitation (a C=0 stretching quantum, or aC — O 
exciton or a vibron) with typical energy of about 1665 cm^' is coupled to lattice phonons with energies 
that are more than one order of magnitude smaller. Due to the assumed strong exciton-phonon coupling, 
the bear exciton becomes self-trapped and forms a polaron, i.e. a new entity that is localized and has 
lower overall energy than the extended exciton. Depending on the relative values of the exciton hopping 
term, the phonon frequency and the exciton-phonon coupling one may arrive to a polaron that is small, 
i.e. localized to essentially on one site, large, i.e. much larger than few sites, as well as having distinct in- 
termediate features depending on the parameters. The Davydov soliton is a special kind of large polaron 
that is formed when phonons respond in an organized, coherent way to the presence of the exciton. Both 
the large polaron and the Davydov soliton are approximate, semiclassical solutions since their creation 
involves a large number of phonons; they are described mathematically through the celebrated Nonlin- 
ear Schrodinger Equation (NLS), an integrable nonlinear partial differential equation f23'2]. Discrete 
breathers, on the other hand, are localized solutions of discrete nonlinear equations and differ substan- 
tially from the extended solitons or large polarons. They involve a local lattice oscillation that is stable 
due to the disparity of its frequency to the linearized frequency modes of the lattice. Although polarons, 
solitons and DB's are in many ways related, especially in some limits, it is useful to differentiate among 
them in order to obtain a clearer understanding of the problem. An important mathematical difference 
between polarons on one hand and solitons or DB's on the other is that the former appear in coupled 
systems involving two fields, e.g. excitons coupled to phonons, whereas the latter are typically single 
field nonlinear equation solutions. In the latter the second field has been eliminated through an approx- 
imate procedure and replaced by an effective nonlinear term acting explicitly on the degree of freedom 
of interest. 



In this review will attempt to describe the various theoretical approaches used in the analysis of the 
acetanilide study and related them to the different experiments. We will start (section 2) with a brief 
description of the original infrared absorption experiment of Careri et al. that linked for the first time 
the anomalous amide-I band to a selftrapped state. In order to proceed with the theoretical analysis 
we will introduce the Holstein Hamiltonian and focus first on a semiclassical treatment (section 3). 
This will enable us to describe the standard adiabatic-like polaron acetanilide picture as provided by 
Scott 110|. Subsequently (section 4) we focus on a fully quantum mechanical approximate analytical 
approach that will give further insight on the acetanilide polaron. In section 5 we discuss the issue of the 
free exciton as appears in the semiclassical as well as the approximate quantum mechanical treatments. 
Subsequently (section 6) we discuss exact numerical solutions of the quantum problem and point out 
similarities as well as differences to the semiclassical solutions. In section 7 we present the pump-probe 
experimental results and compare them to the exact numerics. Finally in section 8 we summarize this 
work and conclude. 
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Fig. 1. Comparison of (a) crystalline acetanilide; dots denote the two hydrogen bond chains formed, and (b) an 
a-helix that is part of a polypeptide. The a-helix forms three hydrogen bond chains. 



2 The acetanilide story 



2.1 Acetanilide and the early experiments 



Crystalline acetanilide (CH3-COONH-C6H5)„ is an organic molecular crystal that forms two close 
hydrogen-bonded chains running along the b direction of the lattice (Fig. 1). The nearly planar amide 
groups have bond distances comparable to those in polypeptides and, as a result, it may be studied in- 
stead of more complex protein and still provide useful polypeptide dynamical information. The infrared 
absorption spectrum shown in Fig. 2 shows two main peaks in the amide-I spectral region; a strongly 
temperature dependent line at approximately 1650 cm^' as well as basically temperature independent 
one at 1665 cm^' II7I8I . While the "normal" 1665 cm^' peak was assigned to one quantum of a C=0 
vibration of an acetanilide hydrogen bonded chain, the origin of the anomalous band was not so clear 
and could be a signature of a Fermi resonance or some structural transition. Both these possibilities 
were ruled out by Careri et al. favoring a then "unconventional" explanation that related the band to 
a "Davydov-like soliton" arising through coupling of the amide exciton to phonon modes |7|. Careri, 
Scott and collaborators used experimental indications and assumed that the C=0 excitation is coupled 
to an optical stretching mode of the hydrogen bond connecting carbon and nitrogen atoms within the 
peptide group |8 9|. Absence of involvement of acoustic modes in the C = O selftrapping was indi- 
cated later through neutron scattering experiments by Barthes et al. 11241 . The theoretical picture that 
the anomalous band is due to a dynamically localized selftrapped state permit the quantitative expla- 
nation of the line temperature dependence using the physics of color centers (25"]. The experimental 
and theoretical picture that emerged in the late eighties asigning the anomalous spectral line to a small 
Holstein-like polaron has been corroborated by the recent experiments of Hamm et al. 1 19 20 21 22] 
that will be reviewed in more detail later. A brief historical exposition of the acetanilide developments 
is given in the Appendix I. 



4 




1640 1660 1680 
Frequency [cm'^] 

Fig. 2. IR spectroscopy of acetanilide. The "normal" and the "anomalous" peak. 



2.2 The Holstein Hamiltonian 



Davydov's molecular soliton model involves the coupling of a vibrational exciton to acoustic phonons. 
The interaction of the excitation with phonons, modifies the latter that in-turn affect the former in a 
selfconsistent way, leading, under several approximations, to an NLS-like soliton 1 3 4 5 23 1 . In the case 
of acetanilide, the initial indications were that the phonons responsible for selftrapping were optical in 
the frequency range 20 cm^^ < CO < 100 cm^'; as a result, Careri et al. used the Holstein Hamiltonian 
instead of the Davydov one. While the theoretical work of Alexander and Krumhansl 1 14] used coupling 
to acoustic phonons, subsequent experimental work by Barthes did not indicate a strong presence of an 
acoustic branch in the low ACN IR spectrum II24I16I17 I . Additionally, recent pump-probe experiments 
seem to show that a unique phonon of frequency approximately equal to co = 50 cm^ ' is responsible for 
the selftrapping phenomenon 1 19|. Although one may not consider the absence of an acoustic branch or 
even an optical branch with dispersion as a closed issue, we may, following the overwhelming number 
of approaches, make the assumption that the physics of the amide-I lines may be described by the 
following Holstein Hamiltonian ll26ll : 



H — Hex + Hph + Hint 

Hex = mf^ (b]bj + 1/2) - y f 



N 

L 



Hph = ll(oY, (b/bj + 1/2 



Htn,=xf^B;Bj(b;+bi) 



(1) 

(2) 
(3) 
(4) 
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with 

where fit ( Bj) creates (annihilates) a C=0 vibron at site j (of sites in total), qj, pj are the position 

and momentum of the y-th optical oscillator respectively while j, bj are the corresponding phonon 
creation and annihilation operators. In Eqs. ([T]|4]i hQ. is the energy of the bare amide-I exciton equal to 
HQ. w 1565 cm^^ while, the optical phonon energy is, as noted previously, h(x) w 50 cm^'. The exciton 
nearest neighbor overlap is of the order / w 5 cm^^ while exciton-phonon coupling is x ~ 25 cm^^. We 
note that acetanilide presents a unique problem in molecular crystals were the three parameters entering 
in the Hamiltonian, viz. J, % and co are known to reasonable confidence. 

The C=0 dipole-dipole interaction transfer integral is very small in acetanilide, and thus, the bear 
exciton band is very narrow and the exciton is immobile. Since the phonon frequency co is much larger 
than J, phonons react very fast to the slow excitons that tunnel from site to site; the phonons then 
follow adiabatically the slowly tunneling exciton motion. Finally, since the exciton-phonon coupling 
% is relatively large compared to J, one expects that a small, very localized polaron may form that is 
not very mobile. In order to treat the ACN parameter regime in the Holstein model we must resort 
in principle to a fully quantum mechanical treatment since we are in the regime where the phonon 
frequency is much larger than the exciton transfer J. Since the first treatment was done through the 
Davydov ansatz, we will first resort to semiclassics and subsequently investigate the connection of this 
approach with the fully quantum mechanical one. 



3 Semiclassical Holstein polaron in one dimension 
3.1 Classical phonon approach 

When either J or % are large many phonons are excited and the zero point quantum motion is not 
important; in these cases we may use the semiclassical limit and treat phonons classically [27|. In many 
practical cases, the excitation of even few phonons is sufficient for rendering the modes classical. Under 
this assumption, while the excitons retain their full quantum nature, the phonon operators b^j and bj 
become c-numbers denoted by b*j and bj, respectively. Furthermore, if in this multiple phonon regime, 
the atomic motion becomes more sluggish, then we may consider the now classical phonons as "slow", 
leading to pj w 0, i.e b* w bj, i.e. bj is real. With these assumptions, the semiclassical Hamiltonian 
reads 

N r- _ 

He,= '£ haB]Bj^J{B]Bj+i+B]Bj^i)+h&b*bf + xB]Bj{b*+bj) . (7) 

7=1 

In order to apply a variational procedure we need first to express the Hamiltonian operator in a specific 
representation, since otherwise operators would have to be equated with c-numbers. Using the one- 
exciton state 

= f V,b]|0)« (8) 
we obtain the expected value of the Hamiltonian of Eq. ^ with respect to the state of Eq. (O: 

Hsc = {MHe.\w) = t [n£l\^l,/-J{^l,*^lfJ+, + +;icofe}Z7; +x|v/(^* +^)] (9) 

Variational minimization of Eq. (|9]l with respect to the b* , gives 



(10) 
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leading, after substitution to the Hamiltonian (|9| to 

(11) 



Hsc — 

.7=1 



hn\^VJ\'~J{YJWJ+l +¥./¥;-i) - £l¥;i' 



Use of Hamilton's equations of motion for the variables \\tj and i\\t*j results in the dynamical equation 

h03 



m\i/j = m\\fj-j{\\fj+i - (12) 



This is the celebrated Discrete Selftrapping Equation (DST) or Discrete Nonlinear Schrodinger Equation 
(DNLS) |28 27|. It is a semiclassical equation for the vibrational exciton of an amide-I mode obtained 
after the complete elimination of the phonon degrees of freedom; it provides an approximate description 
for the dynamics of the exciton when coupled to phonons. 

Within the classical phonon assumption, it is possible to improve the DNLS approximation of Eq. 
( fT2] l by retaining the dynamics of the classical phonon modes. To this effect, instead of performing 
the variational minimization to Eq. (O, we may use Hamilton's equations with respect to phonons and 
obtain a coupled set of two equations, one for the exciton amplitudes and the second for the classical 
phonons |,27J . The resulting set is compatible with the Born-Oppenheimer approximation: 

!^Y^ = -(¥.7+1 +¥;-i) + Cf-^jUj (13) 
^ + «,--a|v,|2 (14) 



where uj a dimensionless displacement with x = cof, a = X^/l/TmJ and y = /ia)/7 [27|. The parameter 
y controls the relative time scales of phonons compared to J. For stationary phonon displacements we 
recover the DNLS equation; this may also accomplished in the limit y ^ 0. 



3.2 Coherent state treatment 

It is instructive to use an alternative approach based on coherent states in order to investigate the semi- 
classical limit of the Holstein Hamiltonian. We now make no assumption on the nature of the phonons 
but assume that the latter are distributed in coherent states, viz. the individual quantum vibrational os- 
cillators are displaced. As a result, we use instead of the state of Eq. dS) the following: 

l¥) = L¥;V|0)«|^-) (15) 

;=i 

\^j)=e^i^r^''j^i\Q) (16) 

where are the complex coefficients that characterize the coherent state at the j-th oscillator and 

|0) is the phonon vacuum. Application of state of Eq. (fTSl l to the Hamiltonian ([T]|4|i, and ignoring the 
constant terms arising from the zero point motions leads: 



N 

Kc = E ^^i^l¥;f --/(¥/ ¥;+l +¥/¥i-i) +^«|^,f +xl¥il'(^i 

i=i 

Variational minimization of Eq. ([TtI i with respect to (|)* gives the condition 



(17) 



that, upon substitution to Eq. ( [TtT i results in the semiclassical Hamiltonian of Eq. ( fTTT ). In other words 
this method leads to the same DNLS Eq. ( fT2] ). as the classical phonon approach. Furthermore, treating 
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(|), /(])*, as conjugate dynamical variables and applying Hamilton's equations to the Hamiltonian of Eq. 
dTTI l results straightforwardly in the equations of motion (I13I14I) . 

We observe that both methods used, viz. the one that treats the phonons directly as classical and 
the other that is fully quantum mechanical but assumes the phonons are in local coherent states, lead 
to identical dynamical equations for the excitons. These are the DNLS equation, when the phonons are 
completely eliminated, or the coupled set of Eqs. (113114b in the Born-Oppenheimer approximation. As 
a result, the variational coherent state method is simply an equivalent but different way of performing 
the semiclassical approximation. 

3.3 Acetanilide analysis 

The semiclassical equations (fT2] l as well as ( 113114b may be used for the analysis of the acetanilide 
problem. We follow Careri et al. |8 | and observe that DNLS equation of (fTSl i has two types of solutions: 

(i) In the limit N °° the extended, free exciton solution = l/^/Ne^'^'^""^^'^^'^ with energy E, 
wavevector k and where a is the lattice spacing. The energy of this solution is 

E,j,{k) ^hQ.-2Jcos{ka) (19) 

where k is the wavevector and a the lattice spacing. The extended solution describes a band of Bloch 
wave excitons with very small dispersion and bandwidth 47; the lowest energy state is that of the \k = Q) 
exciton state, the one that is typically excited in spectroscopy with energy E^x ^ E{k ~0) ~ Ti£l — 27. 

(ii) The selftrapped solution where the probability amplitude is centered at a given lattice site and 
decays exponentially around that site. The selftrapped state may be constructed numerically through an 
iterative procedure starting from the anticontinuous limit at zero coupling, i.e for 7 = 0; it corresponds 
to a simple discrete breather of the DNLS equation|27|. To lowest order the selftrapped state energy is 
calculated by substituting \|/,- = exp{—iEpi,t/h) and ignoring the inter-site coupling 7 that is small; we 
obtain Ept, — TiQ. — x^/h(a. The binding energy then of the selftrapped solution, i.e the total excitonic 
energy gain of the selftrapped solution with respect to the free exciton solution is 

Eh = Eex -Ept~-^+2J (20) 

This is approximately the amount of energy one gains by forming a localized state in one dimension 
compared to the extended, long wavelength Bloch state. 

According the interpretation of Careri et al.|8 10], the temperature dependent peak of acetanilide 
corresponds to the selftrapped state, the peak at 1665 cm^' is the free exciton, while the difference in 
energy between the two peaks gives the selftrapped state binding energy. We note that the selftrapped 
state is a small semiclassical Holstein polaron described trough the DNLS equation; as a result we may 
also call it a discrete breather|27|. As expected, this is a fully localized state that can only acquire 
band character if the eliminated quantum fluctuations are put back in the picture. The semiclassical 
picture for the IR absorption experiments of the amide-I mode assumes that only one optical phonon 
mode of energy hv) « 50 cm^' participates the process; this leads in the adoption of the Holstein model 
for the description of the physics. Three assumptions are subsequently made, viz. (i) the phenomena 
may be described by a one dimensional model due to the quasi-one dimensional nature of the crystal 
in the direction of the hydrogen bonds, (ii) multiple phonons are excited, and as a result, the phonon 
variables may be taken to be classical and (iii) these classical phonon variables may be very slow. The 
semiclassical picture then gives a clear physical picture for the two peaks assigning the "anomalous" 
one to a small polaron selftrapped state and the normal one to a free exciton state. 

Although the semiclassical picture is successful in the assignment of the two spectral lines it cannot 
address the temperature dependence of the anomalous peak; for the latter one needs to consider directly 
phonon excitations, that nevertheless have been eliminated due to the classical phonon approximation. 
One way to proceed is by ignoring polaron hoping motion altogether (since 7 is small) and requantizing 
the previously assumed classical phonons. We obtain the picture of Fig. [3] where the two parabolas 
correspond to phonons around the no-exciton and one-exciton states respectively; the latter state is 
displaced with respect to the no-exciton state due to the exciton-phonon interaction. In this picture, the 
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Phonon Coordinate q 

Fig. 3. The displaced oscillator picture; Epi, is the exciton energy gain due to coupling to phonons. 



polaron line temperature dependence is obtained similarly as in the cases of color centers in molecular 
solids [25 1. It is provided through the sum of all overlaps between equal phonon quanta states of the 
ground and displaced excited states weighted by the initial (thermal) population of the phonon states, 

viz. m 



W{T)=Y^P„\{n\^„) 

;i=0 



exp(-n — 



(21) 
(22) 



where \n) is ground state with n phonons, |(|)„) is a displaced oscillator state with equal number of 
phonons, while P„ gives the equilibrium distribution of phonons at temperature T . The Frank-Condon 
factor W{T) is found to be [lOJ: 



W{T) = exp 



coth 



/ico 

7k^ 



i-^csch — -— 
rico2 V2kbT 



(23) 



where /q is the modified Bessel function of first kind and order zero. It is found that this expression 
describes correctly the temperature dependence of the anomalous peak lITOl . 

Although the semiclassical picture with the quantum mechanical adjustment for the displaced os- 
cillator Franck-Condon lines appears to describe fully the ACN data, alternative explanations and other 
issues have been raised in the past. These include (i) A Fermi resonance, topological defects or non- 
degeneracy in the hydrogen-bond may lead equally to a similar double peak feature in the absorption 
spectrum, (ii) The overall validity of the approximations involving Davydov-like adiabatic derivations 
may not be justifiable, (iii) Can acetanilide be described faithfully through a one dimensional Hamilto- 
nians? In what regards alternative explanations for the double peak, there is consensus that the exper- 
imental and theoretical work performed appears to rule most of them out. The experiments of Barthes 
et al. have ruled out convincingly the topological defect suggestion while also shown that the hydro- 
gen bond non-degeneracy suggested by Austin and coworkers is not valid ll29..30j . Furthermore, the 
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analysis of Careri and Scott as well as the recent analysis of Hamm have shown that the occurrence 
of a Fermi resonance does not seem very likely f8"20 |. The validity of the adiabatic approach used by 
Scott and collaborators has been also discussed, however works focused more on the general Davydov 
soli ton problem rather that the specific acetanilide application II12I31I . Acetanilide is a narrow band 
solid and, as a result, the occurrence or a Davydov-like extended object is unlikely. As a result the ac- 
etanilide small polaron regime may be better analyzed either perturbatively using 7 as a small parameter 
or through the direct numerical diagonalization. Finally the issue of the true dimensionality of the ACN 
selftrapping phenomenon demands also attention since small polarons in one and three dimensions have 
quite different properties ll27l32l . These issues will be addressed in the following sections. 



4 Quantum analysis in the small 7 limit 

While the semiclassical analysis is useful as well as intuitive, the mere fact that the hoping matrix 
element J is quite small in acetanilide allows for a complete quantum mechanical solution in the limit 
7 = and a subsequent use of 7 as a perturbation parameter 1 14 15 1. This approach is fully quantum 
mechanical and thus independent of any type of approximations. 



4.1 Quantum 7 = limit 



When the nearest neighbor coupling is zero the problem can be diagonalized exactly. The Hamiltonian 
becomes for 7 = 0: 

N N N 

H = h^Y. [B]B j + 1 /2) + /zco £ (b]b j + 1 /2) + X £ B]B j [b] + b^) (24) 

7=1 7=1 7=1 

We perform the Lang-Firsov transformation lf33l by introducing the unitary operator U = e^^ where 



S^^yB]Bj(b]-bj 
For the diagonalization it is useful to consider the following identities: 



u'^b]BiU = b]Bi 
U^b]u = b]~ ^BiBi 



U^bi^'biU = b]bi - ^B]B,{bJ+bi) 



/iCO 



\Tl(£>J 



U%iU = bi- i^B]Bi 



b;bib;b, 



The transformed Hamiltonian H = U^HU after removing the zero-point motion terms becomes 



N 
7=1 



no. 



h03 



B/Bj-^B/B/BjBj + hcob/bj 



(25) 

(26) 
(27) 

(28) 
(29) 

(30) 



Equation ( [30] l may be diagonalized directly, leading in the case of one vibron to the energy spectrum: 

E=(nQ.-—\+nn(a (31) 

for « = 0,1,2,.... We note that the one vibron excitation has lowered its energy due to the coupling to 
the phonons by an amount equal to (see Fig. [3]l 



h(£)' 



(32) 



The quantity Eph is the polaron binding energy determined quantum mechanically in the 7 = limit. 
This result coincides with the semiclassical polaron binding energy calculated in section (|3.3l l. 
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4.2 Quantum J^O limit 

When J y^Q the unitary transformation appHed to the vibron transfer term gives 

N . , 

U'^HjU = JY. [U^B]Bj+iU + u''B]Bj-iU) = 



^ / 



N 
I 

7,,+,=7exp{-^ 



+ iBjBj+i +Jjj^iBjBj^i 



(33) 
(34) 



In order to proceed we assume first that the system is at zero temperature and the phonons at each site 
are in their corresponding ground state, viz. |0) = |0)p/, = ...|0)j-i |0);|0)j+i...; using this state we find 
an effective coupHng between adjacent sites: 



= (0|-/;,;+l|0) =yexp 



\hoiJ 



(35) 



The effective bandwidth of the polaron tunnehng is thus reduced as a resuh of the interaction with 
phonons |34|. This renormahzation of the polaron bandwidth is an approximate, mean field result, that 
ignores local phonon fluctuations and has a small effect when J is small (J ^ hd)). While the Lang- 
Firsov transformation does not depend on the value of J, the averaging over the zero phonon states 
should be done for small hopping matrix elements. 

Using Eqs. (l35T l we may now write the transformed Hamiltonian in the 7^0 limit for the case of 
zero phonons as follows: 



Bi^Bi-J, 



{b/Bj 



+i+Bj'Bj^i)-^Bj^B/BjBj 



(36) 



If we assume that only one C — O vibron is present while phonons are in their ground state, i.e. 



l¥) = LV;B/|0>«|0)p/, 



(37) 



then we can directly diagonahzed the tight-binding Hamiltonian. Using Eq. ( |37] | in Eq. ( l36b leads to 



The energy spectrum for the one-vibron sector is 



E{k)^[ no. 



TlK) 



-2Je~^'^' cos{ka) 



(38) 



(39) 



The approximate small-7 procedure leads to a small polaron band that has gained energy ~x^/h(£i with 
respect to the center of the bear exciton band (obtained for X = 0) and has a reduced bandwidth equal 
to 47e^'ra) . This polaron band is the only one-exciton, zero-phonon solution in the 7^0 limit. The 
approximate energy £ (A;) for the small polaron band of Eq. (|39l ) is identical to the one obtained by Scott 
using degenerate perturbation theory ifTOll . 

Let us briefly summarize the findings of the approximate but fully quantum mechanical treatment: In 
the search for the eigenstates of the Holstein Hamiltonian, the exact quantum mechanical calculation at 
7 = results in an eigenstate that has energy lower by an amount equal to Ebp — —X^/hu) from the bear 
exciton state. This state is thus favored energetically over the bear exciton that exists when X = 0. The 
value of the binding energy of this new polaron state coincides with the semiclassical (DNLS) binding 
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energy. For exciton hoping J small we find (at T = 0) that this state acquires dispersion with bandwidth 
B = 47exp [— (x/Zzco)^] . Consequently, this solution provides the small polaron band formed as result 
of exciton-phonon coupling and the small J value. Since the semiclassical and quantum mechanical 
polaron binding energies are identical, we may consider that the localized semiclassical polarons begin 
to tunnel when quantum fluctuations are included to low order and form a band with bandwidth B. We 
point out that for acetanilide B w 16exp(— 1/4) cm^^ « 12.4 cm^', i.e. the bandwidth reduction due 
to exciton-phonon coupling is small. 

Although the polaron solutions in the semiclassical and quantum treatments coincide, the free exci- 
ton solution found in the context of DNLS seems to be absent in the quantum mechanical analysis. It is 
in fact easy to see that in the fully quantum approach the exact free exciton state is not an eigensolution 
in any dimension since any extended exciton wave state generates phonons through the exciton-phonon 
coupling term in the Hamiltonian. Since the presence of the free exciton peak is central to the explana- 
tion of the acetanilide amide I spectrum one must find a way to bypass this problem. One possibility 
that is compatible with the semiclassical picture is to consider that the free exciton state is actually an 
excited polaron state that involves a number of excited phonons. However since the phonon involved in 
the process has energy equal to 50cm^^ while the polaron binding energy is less than half of this value 
(approximately 1 6cmr ' ), an excited polaron plus one phonon cannot match the value of 1 6cmr ' , except 
for exceptionally large values of % that are not reasonable for molecular crystals. Another possibility is 
that the description of acetanilide as a one dimensional solid is not sufficient. It is known semiclassically 
that in 3Z) the polaron state may be separated by a barrier from the free exciton state and, as a result, 
both states may coexist in some fashion [27J . If such a picture survives in the fully quantum case and 
for the proper parameter regime, it might then be possible to identify the two states accordingly. These 
issues will be addressed in detail later. We note here that the small-/ analysis can be done easily in three 
dimensions leading to a small polaron band with energy 

^(k) = (^Q. - - 2e"('^'^ [/,cos(Mx) +-/vCos(/fcyflj,) +7,cos(/fc,fl,)] (40) 

where J = {Jjc,Jy,Jz) are the hopping rates for the three Cartesian directions respectively, k = {kx,ky,k^) 
is the wavevector and a = {a;t,ay,a^) the lattice spacings. 



4.3 Brief summary on various approaches 

It is worthwhile to summarize briefly the various equations obtained so far through three different ap- 
proaches. 

- The semiclassical DNLS equation may be obtained variationally either by (a) considering phonons 
classical or (b) assuming that phonons are quantum but distributed according to coherent states. In 
both cases we obtain identical results, viz. energy gain for the polaron and no bandwidth reduction. 
In the Born-Oppenheimer approximation both approaches lead to identical dynamical equations as 
well. 

- When we solve the Holstein model fully quantum mechanically but in the small hopping limit we 
find a polaron band with energy gain and a reduced bandwidth. 

Comparing these approached we find that all three of them (pure semiclassical, phonons in coherent 
states as well as quantum mechanical) give the same polaron energy gain, yet, the bandwidth reduction 
does arise only from the fully quantum mechanical approach. Furthermore, the free exciton state is 
a stationary solution only semiclassically. While the semiclassical approaches result to an effective 
nonlinear equation of motion, viz. DNLS, the approximate quantum methodology leads to a QDNLS 
Hamiltonian 1351 . The latter provides a fully quantum mechanical reduced description for the excitons 
provided the hopping rate J is small and phonons are not excited on average from their ground states. 

From Eq. (l36T l we see that the on site term of the QDNLS Hamiltonian contains already the polaron 
binding energy shift while the hopping term is also renormalized. When only one exciton is present, this 
Hamiltonian describes a small polaron band. When more exciton quanta are present, it may describe 
multiexciton states as well as interaction among them. It is clear that the QDNLS approach is valid 
when intersite phonon correlations are not important. 
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5 The free exciton state in the semiclassical analysis 

In order to focus deeper into the double peak acetanilide structure we may investigate in more detail the 
semiclassical picture before ultimately probe directly the fully quantum regime using exact numerics. 
In this section we will look into the free exciton state both in the Born-Oppenheimer approximation in 
one dimension but also treating the three dimensional case. 

5.1 Born-Oppenheimer approximation 

In the fully adiabatic picture the classical Einstein oscillators of the Holstein Hamiltonian have no dy- 
namics but follow the exciton motion. We may still consider these oscillators classical but allow them a 
characteristic time scale for response to the exciton motion. In this case we need to solve the set of equa- 
tions ( 113114b that determine the stability of the adiabatic solutions. The acetanilide values correspond to 
a w 2.3 and y w 10; thus if we substitute the semiclassical adiabatic solutions for the small polaron and 
free exciton of section (2.3) and use Eqns. (113114b we will find what is the regime of stability of these 
solutions. The adiabatic polaron state may be constructed via a numerically exact procedure from the 
anticontinuous limit; its linear stability has already been discussed f27 |. The stability of the free exciton 
state, on the other hand, may be assessed numerically by direct substitution of the state in the equations 
of motion of Eqs. ( I13I14I) ; the result is shown in Fig. (0]). In order to engage the vibrational lattice to 
the free exciton we have added initially a small component of uniform noise in the classical phonons. In 
the figure we show the free exciton evolution in time units tex =Y^/5 = Jt /5; we observe that after few 
exciton oscillations the free exciton becomes unstable and forms localized structures. The phenomenon 
of the semiclassical polaron formation from the free exciton state is controlled by a characteristic po- 
laron formation time Zp. This waiting time for the onset of localization depends strongly on the specific 
parameter regime as well as initial amplitude to the phonon excitation. Extensive simulations show that 
Tp — > oo when Y ^ 0, i.e. when we address the fully adiabatic regime. In the opposite limit of y > 1, on 
the other hand, the semiclassical free exciton becomes very quickly unstable and localizes. We expect 
that as the temperature of the system increases, this effect becomes more dominant leading to a collapse 
of the free exciton into localized states. 

The use of a finite parameter y — h(o/J in the semiclassical problem introduces a new time scale 
in the Holstein Hamiltonian; while T is the natural time scale for the phonon evolution (measured in 
phonon periods) , the exciton time is measured in Zex = T^- In the adiabatic regime y <C 1 and in the 
course of many exciton oscillations there are very few phonon circles. Using in Eqs. (113114b as initial 
condition a free exciton state, we observe that the slow phonon change simply modifies the local site 
energies of the exciton equation; as a result the free exciton suffers very weak phonon perturbation and 
survives for long times. In the opposite regime where y^ 1, many phonon oscillations occur in few 
exciton ones; as a result the interaction term acts as an effective multiplicative noise term (for some 
initial phonon conditions) leading to a very rapid collapse of the free exciton into localized states. This 
purely classical effect is enhanced in the fully quantum regime through the zero point motion. 

Since the parameter regime of acetanilide is in the range y > 1 we find that the free exciton is 
generally unstable to even quite small perturbations. This behavior poses a problem to the standard one 
dimensional adiabatic explanation since, for this approach to be valid the free exciton state should be 
stable as is the small polaron. One way to proceed within the semiclassical context is to consider the 
possibility that the selftrapping phenomenon takes place in three dimensions. Since, as noted earlier, in 
3D a barrier separates free exciton and polaron, it may be that the barrier helps in stabilizing both states. 
||27I . We will proceed with the semiclassical analysis of the 3D Holstein model in the next subsection. 

5.2 Semiclassical 3D problem 

It was pointed out recently that the dipole-dipole interaction responsible for the nearest-neighbor transfer 
of excitons in ACN extends not only along the hydrogen-bonded chain direction but also in the plane 
perpendicular to it l[36l ; this point will be detailed in sections (|6]l and ( 17.1b . In order to perform the 
semiclassical analysis for the general problem in three dimensions we write the semiclassical Holstein 
Hamiltonian as follows: 
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e=< lattice site 



Fig. 4. Time evolution of the site probability as a function of exciton time tex = t/Sy = Jt/5 for 7= 10 and 
a = 2.3. We used a periodic one dimensional lattice of 100 sites. The initial condition is that of a free exciton 
of zero momentum and uniform noise in the initial (dimensionless) phonon displacements ufs with amplitude of 
the order of 10"^ (the phonon lattice evolution is not shown). The free exciton is initially stable but subsequently 
collapses to localized states. 



= £ [^iniVjl' - (Av)j + to^l^j +Xl¥jl'(^| ] (41) 
j 

(43) 

where j = {jxjyjz) h^JyJz the transfer rates in the three perpendicular axis directions. We 
assume that the hydrogen bonded axis coincides with the x-axis ; the parameters Jx, Jy, Jz may take pos- 
itive (negative) values leading to negative (positive) hopping rates. We may now perform the variational 
minimization with respect to the phonon c-numbers b*^, bj and obtain the variational energy Ey. 



r, 



-Jx{WM+l,jyj,-^j\^-2) -^(|V;.J,+l,A-¥jP-2) -Jz{\^Mjy,h+l -¥jl^-2)] 



(44) 



where d is the dimensionality. In order to probe into the variational energy E^, we make the ansatz 

=A^I-'^lril-'>lel-''^l (45) 

where 1^1 < 1, |ri| < 1, |0| < 1. The three parameters ^,ri, and 9 take positive (negative) values when the 
hopping overlap in the corresponding direction is negative (positive). The normalization factor is found 
to be 



2\ 1/2 



1 



2\ 1/2 



1-9^ 
1+92 



2\ 1/2 



(46) 
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Substitution of the trial function of Eq. (05]) into Eq. (l44l i we obtain the following expression for the 
variational energy 



E, = -4 



-J. 



11 



l+ri2 



e 



1+92 



(1 - + (1 + V) (1 -e2)(i + 04) 



/iCO 



2\3 



(1+02)3 



(47) 



In the fully symmetric case of Jx — Jy — = J > v/e have — r\ — Q = r\ and obtain the general 
expression d dimensions ll27l 
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^2 (i_ri2)^(l+riY 



l+ri2 h(0 {l+y\^yd 



(48) 



Figure |5] shows the variational energy (in units of 7) in ID and 3D for the isotropic case with all 
negative hopping overlaps (J > 0). In ID, the curve is barrierless for any set of parameters, regardless 
how small the coupling is. That is, a partially localized solution r| < 1 has always smaller energy than the 
extended solution r| = 1, and furthermore, this solution can be reached without having to surmount any 
barrier. The behavior in 3D is drastically different. For couplings below the critical threshold value Uc = 
y^lx^./hoxf = 3.288, the extended solution with r| = 1 is the lowest energy solution. For couplings larger 
than this critical value , a strongly localized state r| ^ 1 becomes the lowest energy solution 1271 . 



5.3 Polaron and free exciton coexistence in three dimensions 

The effective semiclassical barrier between the free exciton and small polaron forms when a = « 
3.288. When a < the minimum energy state is the free exciton while for a > ttc is the small polaron. 
When the parameter values are such that an adiabatic double well structure is formed, the transition 
from free exciton to small polaron occurs through quantum lattice fluctuations. A simple way to include 
the latter while retaining the intuitive adiabatic picture is by treating r| as an effective "configurational 
variable" and quantizing the double well potential problem. In other words, we may consider that the 
dynamics of the transition between polaron and exciton may be described by a particle of effective mass 
m* that is moving in the adiabatic potential. Upon diagonalization of this effective one dimensional 
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Fig. 6. (a) The adiabatic potential for the 3D polaron as a function of displacement; the horizontal lines correspond 
to the first two energy eigenvalues while dots the next two for a = 4 and m* = 5. In (b) the spatial distribution of 
the probability for the lowest three eigenfunctions, the ground state (GS), first (1) and second (2) excited states, 
respectively. 

problem we may obtain the states involved and understand their dynamics. The results are shown in Fig. 
(ID) for a = 4 where we portray the adiabatic potential, the four lowest energy eigenvalues as well as 
the three lowest eigenfunctions. The equilibrium position is zero for the undisplaced free exciton state. 
For the parameter regime shown the lowest potential minimum corresponds to the small polaron while 
the smooth, wide secondary minimum is the free exciton minimum. 

From the numerical solution of the spectral problem we find that the two lowest states are relatively 
close in energy and form a type of doublet. The ground state wavefunction has a single peak that appears 
to have two components. While it is centered close to the polaron minimum it has a considerable overlap 
with the exciton minimum as well. The first excited state, on the other hand, is doubly-humped and 
centered in the wide, shallow secondary minimum of the free exciton but with significant amplitude in 
the polaron minimum. The ground state-excited state doublet has distinct polaron-like and exciton-like 
structure respectively. If the 3D Holstein model is a good representation of the acetanilide problem, then 
it is possible that these are the two states that are relevant for the explanation of the double peak in the 
ACN absorption spectrum. 

In the previous analysis we used a relatively large value of the nonlinearity parameter so that we 
are withing the double well regime of the adiabatic potential. This choice may be justified in the fol- 
lowing way: Since we know that in reality the polaron hopping rate is reduced we may in an ad-hoc 
fashion use the effective rate of Eq. (l35T l instead of just J and obtain a new nonlinearity parameter 
deff = exp [5(7^)^] OC that is larger than a. This way the effective nonlinearity increases for the same 
acetanilide values and the system is found to be closer to the adiabatic double well regime. Clearly 
further investigation on this approach is necessary. 



6 Numerically Exact Diagonalization of the Holstein Hamiltonian 

With recent advances in computer technology it became possible to diagonalize the Holstein Hamilto- 
nian on a numerically exact level, even in the 3D case li37i38J . This allows one to study to what extent 
the conclusions from the DNLS (e.g. the existence of a barrier in the 3D case) manifest itself in the full 
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quantum regime. To that end, a concept introduced by Trugman and coworkers is used II37I38I . In brief, 
Bloch states with momentum k are constructed: 

I' = ^ y <i>e"'^ (49) 

with site states <I> that are identical for each site j. The site states are expanded in a basis cp{„;}: 

* = E '■{",}^{«,} (50) 
{",} 

with (in one dimension): 

9{n,} = |lo;---,n-l,no,n+l,---) (51) 

where the ' Iq' denotes that we have one quantum in the exciton coordinate at site (since the Hamilton 
is quantum-conserving with respect to the exciton coordinate, it is block-diagonal and we consider here 
only the one-quantum manifold of states). The «, are the numbers of phonon excitations at site / around 
the exciton at site 0. 

The essential trick of Refs. f37'381 is an efficient strategy to select basis states, which we review here 
only very briefly. We start from a 'root' basis state without any phonon excitation: |lo; ...,0,0,0, ...). For 
each generation (total number of generations: Nh), basis states are added corresponding to jumps along 
one of the off-diagonal elements of the Hamiltonian Eq. (|2]B1) in a tree-like manner. This can either be 
one phonon jump with a matrix element 

(lo;...,n-i,no,n+i,---|^/^|lo;---,n-i,no + l,n+i,---) = - ^"0 + IX (52) 

or a shift along one of the crystal directions with matrix elements 

(lo;.. no, n+i, •••1^1 lo;---,n-2,n-i, no, •••) =Jexp{ika). (53) 

In optical spectroscopy, where a A: = selection rule exists, the phase factor exp{ika) disappears. The 
basis gets a pyramidal shape with maximum Nh phonon excitations at site to one phonon excitation at 
one site ±{Ni, — 1) in either lattice direction. The basis is infinite with respect to lattice transformations 
and fully accounts for the translational symmetry of the crystal (Eq. |49] |. However, the basis is finite 
with respect to the phonon-exciton distance, and the phonons are confined to the exciton by construct 
within ±(A^/, — 1). The number of basis states scales as (9 ((Z) + 1)'^''), where D is the dimensionality of 
the system. In 3 dimensions with generation number A^/, = 11, this reveals as many as 1.4 x 10^ basis 
states. Nevertheless, the matrix is sparse and can be partially diagonalized using a Lanczos algorithm. 
On an AMD Opteron processor with 12 GB main memory, the ground state is calculated within minutes, 
or, if the whole absorption spectrum is needed, within typically a day. 

Before we start, one comment is in order: In order to distinguish a free-exciton from a polaron 
solution, we will use the phonon displacement q asa criterion. It might seem more logical to investigate 
the extension of a solution in real space, similar to in the DNLS case (Sec. l5.2l i. However, even though 
the free-exciton would, of course, be perfectly delocalized in real space, any polaron solution would be 
as well in the full-quantum case and one could not distinguish between polaron and exciton in this way. 
This is by construct of the Bloch ansatz Eq. ( |49] l. which accounts for the translational symmetry of the 
problem. In the strict meaning of the word, there is no true, spatial self-trapping in the full-quantum case. 
This is fundamentally different from the semi-classical DNLS case, which does break the translational 
symmetry of the problem and reveals solutions that are localized in real space. Only in the limit of 
infinite effective polaron mass, where the polaron-dispersion relation would becomes infinitely flat, one 
could construct spatially localized wavepackets, which then would indeed be eigenstates of the polaron 
Hamiltonian. Comparison with Eq. ( [39] l shows that a flat dispersion relation is obtained for x/Zzco ^ oo. 

Fig. [T] shows the results of a numerically exact matrix diagonalization of the Holstein Hamiltonian 
for exciton-phonon coupling x/hoi = 3.5 with the exciton coupling J/hoi varied from to 3. Shown is 
the projection of the corresponding ground state wavefunction onto the phonon coordinate q. With this 
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Fig. 7. The Holstein Hamiltonian in 3D with isotropic negative exciton coupling: Projection of the ground state 
wavefunction onto the phonon-coordinate q. The exciton-phonon coupling was set to x/^K* = 3.5 in this calculation, 
and the exciton coupling J/hai was varied from to 3. 



set of parameters, the system should be close to the DNLS regime (x/ho) ^ 1), so one expects that one 
can directly compare results from Sec. l5.2l with that from the exact diagonaUzation. 

Indeed, for small exciton phonon couplings J, which corresponds to a large a = 2%^ / hoxf , the 
wavefunction is displaced with respect to the phonon coordinate, close to the value \/2x/li(i) expected 
from the displaced oscillator picture (Sec. 14. 11 1. Referring to Fig.|5] this corresponds to a polaron solution 
with small r| for a > a^, which in this parameter range constitutes the lowest energy solution. However, 
as the exciton coupling is increased to ///ico « 2.2, which for the given exciton-phonon coupling relates 
to the critical value = 3.288, the ground state wave function jumps almost suddenly back to q ^ 
0. From Fig. |5] we expect for a < that the extended state with r| = 1 becomes the lowest energy 
solution, and accordingly, the phonon-displacement q vanishes at this point. If the system were truly 
classical, the switching between the two regimes would occur instantaneously, but since the quantum 
mechanical wavefunction delocalizes, a relatively small parameter regime exists where polaron and 
exciton solutions coexist (similar to Fig.|6]l. However, we note that even though exciton-like solutions 
are found in the quantum case, they are not exactly free excitons (the displacement q is not exactly zero). 
The free exciton is not an eigenstate of the Holstein Hamiltonian. 

Fig. [8^ puts these results into a more general perspective. Shown is the expectation value of the 
phonon displacement (q), relative to that expected from the 7 = case q,„a.\ — V^JC/^K), as a function of 
the two parameters of the theory x/h(0 and J/h(0. The thick line separates the two regimes of the semi- 
classical DNLS, defined by a^- = 2%^. / h&J ~ 3.288: Left from that fine, the DNLS predicts polaron 
solutions with displacements close to the maximal value q„,ax = V2x/h(a, whereas extended states 
with r| = are expected on the right side of this phase diagram. For larger exciton-phonon coupling, 
X/Zico >« 2 the quantum solution essentially follows that line, except that the turn over between polaron 
and exciton solution is not discontinuous (as discussed above). Nevertheless, as x/h(0 is increased from 
2 to 3.5 (Fig.[8}3 middle and top), the transition occurs more and more abruptly. 

In contrast, for x/^ift) <~ 1 the quantum solution deviates strongly from the DNLS prediction. In 
particular, polaron solutions with large displacements (q) may exist even for relatively large exciton 
coupling J, where the DNLS would strictly reveal an exciton solution (Fig.[8}3, bottom). This difference 
between semi-classical and quantum result can be rationalized by taking into account the second term 
ofEq. 



2/exp [-x^/(/!C0)^] cos(fea) 



(54) 
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Fig. 8. The Holstein Hamiltonian in 3D with isotropic negative exciton coupUng: (a) Displacement of the phonon 
coordinate q (in units of x/h(s)) of the ground state according to a full-quantum diagonalization of the Holstein 
Hamiltonian. Results are plotted as a function of exciton coupling J and exciton-phonon coupling both in units 
of the phonon frequency ha. The thick line separates the two regimes predicted from the semi-classical DNLS: 
ttc = \JT^JVmJ — 3.288 (b) Horizontal cuts through (a), exemplifying the quantum case (bottom, x/'^fJ = 0.5), 
the close to classical case (top, x//ico = 3.5), and an turn-over regime (middle, x/Sco = 2.0). The grey lines represent 
the results from an exact diagonalization, the black lines that from the DNLS, interpreting the interaction term 
in Eq. i ll U as an effective coupling Xe// = XlVpJ hence the maximum displacement found in the peak of the 
wavefunction is ^ = \/2xlVmflAp/''^- 1" either case, the displacement q is plotted in units of \/7.x/h(£i. 



which may be considered a quantum-correction factor relative to the polaron binding energy in the 
semiclassical DNLS case (Eq. [STT i. The origin of the quantum-correction factor is exciton tunneling 
that occurs when a spatially localized solution (Sec l4.2l i starts to delocalize due to a non-zero exciton 
coupling 7 ^ 0. In that case, a phonon displacement follows the exciton as the latter hops from one site 
to the next. In fact, the factor exp(— x^/(/iCo)^) equals the overlap integral (Franck-Condon factor) for 
this hopping with the phonon coordinate in its quantum mechanical ground state (see Eq.[35]l. For large 
exciton-phonon coupling, we have exp(— x^/ (^w)^) 0, and Eq. ( [39] l would equate Eq. (ISTl i. In that 
hmit, the polaron dispersion relation would become flat (i.e. the polaron energy would be independent 
on wavevector k), and truly spatially localized solutions, as we obtain them from the DNLS, could 
indeed be eigensolutions of the quantum Hamiltonian. This is the classical limit, where tunneling is 
not possible and the hopping probability becomes zero! However, for small exciton-phonon coupling, 
the polaron energy is lowered by an additional term 2J / [hu))^) , which may render the polaron 
more stable than the exciton, even when the DNLS would predict the opposite. 

It is instructive to compare the results of this chapter with those of the coherent state treatment. We 
have seen in Sec. I3.2l that the product ansatz Eq. [T5]and [T6l directly leads to the semiclassical DNLS 
without invoking any further approximation. The DNLS, however, fails in certain parameter regimes 
where x/h(i> <« 1 . The essential difference between coherent state treatment and the approach discussed 
above is that the latter uses a much larger basis (Eqs.|49]and ISTl ). In particular, the 'basis' of the coherent 
state treatment is local in the sense that the polaron is confined to the exciton by construct. In contrast, 
the basis Eqs. l49land |5T| allows phonon excitations away from the exciton with, e.g. a phonon sitting 
at site 1, while the exciton is sitting at site at the same moment (i.e. «+i = 1). Counterintuitively, that 
phonon at site 1 does not interact with the part of the exciton that is also residing on site 1 when shifting 
<I> (Eq. ISOl l by one site along the lines of the Bloch ansatz Eq. |49] There is no classical counterpart to 
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this phenomenon; it reflects the non-local character of quantum mechanics. In the hmit %lh{i> 3> 1, the 
nonlocal character become less important, and the DNLS becomes a good approximation. 



7 Pump-probe experiments on ACN 
7.1 Experimental Results 

With these introductory words in mind, we now turn to our recent femtosecond pump-probe experiments 
on the C=0 band of crystalline ACN [20 21 1. (The pump-probe response of the NH band of ACN has 
been studies as well 1 191221 . which however shall not be discussed here.) In such a pump-probe exper- 
iment, the various modes are first excited with an intense short IR laser pump pulse, and the response 
upon that excitation is then probed by a second short, but much weaker probe-pulse. The informa- 
tion gain of femtosecond pump-probe experiments, as compared to simple IR absorption spectroscopy 
(Fig.|2]i, is manifold: 

- The anharmonicities of the individual modes can be measured, as sketched in Fig.|9j5: When pump- 
ing a state, its first excited state is populated, and the probe pulse is then probing both the 1-2 
upward transition and the 1-0 downward transition. If the state is harmonic, both upward and down- 
ward transition appear at exactly the same frequency, but with opposite signs, and cancel each other 
exactly. In contrast, if the state is anharmonic, they lead to two bands with alternating signs separated 
by the anharmonicity of the transition. 

- If two states are coupled, excitation of the one state will lead to an response (frequency shift) of 
the other. It shows to what extend the various transitions are entangled, or whether they can be 
viewed as separated problems. If, for example, the two lines in ACN would originate from molecules 
in different conformations or surroundings (topological defects), as has been suggested various 
times 1 39 29 1, pumping of the one state would not cause any response of the other state (in con- 
trast to what is seen in the experiment). On the other hand, in the case of a Fermi-resonance, another 
often discussed alternative explanation of the doublet in the ACN spectrum |40|, a very characteris- 
tic coupling pattern is expected. Investigating such coupling patterns in detail, compelling evidence 
against both of these possibilities could be provided (21]. 

- The lifetime of the pumped states can be measured with femtosecond time resolution. In this way, 
the original idea of Davydov [4J can be tested, namely to what extent vibrational self-trapping may 
stabilize the excitation for timescales that would allow to make the energy available for subsequent 
biological processes. 

Fig.|9h shows the result of pump-probe experiments on the two bands of crystalline ACN fSOl : When 
resonantly pumping the "anomalous" band ( 1 650 cm^ ' ) of ACN, the band bleaches (negative response) 
and a positive band emerges at 1644 cm^^ In contrast, when resonantly pumping the "normal" band 
(1666 cm^'), hardly any bleach of the band itself is observed. Nevertheless, the anomalous band does 
responds with a signal which is similar in shape, but is slightly smaller than when pumping it directly. 

In a somewhat heuristic approach, Edler et al. have discussed the distinctly different responses of 
the two bands of ACN in Ref. |20|. The self-localized state was treated along the lines of Scott's the- 
ory 141], which would relate the anharmonicity of the anomalous band to the —j^/hmBj^Bj^BjBj term 
in Eq. ( l36b . At the same time, it was assumed that one can 'cut out' the exciton part of the Hamiltonian 
Hex (Eq.|2]i from the remainder by writing: 

Hex = Ti^Y. [Bj^^j + 1/2) - / L [B/Bj+i +B/Bj^ij (55) 

The discussion of the exciton is a little more involved: It turns out that, apart from the artificial separa- 
tion from the phonon/polaron part, this Hamiltonian is not a very realistic description of molecules. A 
chemical bond is not harmonic in reality, but is better described by a Morse-potential. Vibrations of real- 
worlds molecules are intrinsically anharmonic (on-site anharmonicity). For example, the anharmonicity 
of a C=0 vibrator, defined as the difference in frequency of the 0-1 and the 1-2 transition, lies between 
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Fig. 9. (a) Linear absorption spectrum (top) and pump probe spectra of the C=0 mode of crystalline ACN at 93 K 
for two different narrow band pump pulses chosen to be resonant with each of the absorption bands. The arrows 
mark the center frequency of the pump pulse, which was spectrally narrow enough to pump only either one of the 
transitions, (b) Level scheme of the system, explaining the distinctively different response of both modes. Adapted 
from Ref. llol. 



10-20 cm . Hence, we should use instead for the exciton Hamiltonian: 

where A is the on-site anharmonicity. Note that the additional term does not affect the one-exciton states 
nor any polaron solution in the one-exciton manifold where B B B jB j = 0, which is why we normally 
can ignore it. It also doesn't affect the outcome of normal absorption spectroscopy along the lines of 
Fig.|2] The term however does affect the way how one potentially sees one-exciton states in pump-probe 
spectroscopy, since the two-exciton manifold is reached in such experiments. 

In light of this discussion, the experimental result shown in Fig. |9^, bottom, is in fact surprising. 
In this experiment, a vibrational state (the higher frequency band) is pumped, but cannot be depleted. 
Similar pump-probe experiments have been done by Hamm and coworkers on hundreds of molecules, 
but a result of this sort has never been obtained. 'Normal' vibrational states are intrinsically anharmonic, 
and as a result, one expect to see a pump-probe response similar to that of the 'anomalous' band, i.e. a 
bleach/stimulated emission at its original 0-1 position and a positive 1-2 absorption frequency shifted 
by the intrinsic anharmonicity A. In Refs. II20I211 . several alternative explanations have been excluded, 
and the only possible explanation left over was: The higher frequency peak in the absorption spectrum 
of ACN behaves as if it were harmonic, despite the fact that it is composed of C=0 vibrators that are 
intrinsically anharmonic. Only if the state is effectively harmonic, upward and downward transitions in 
Fig.lHJ? could cancel exactly and the zero-response seen in the experiment is expected. 

Interestingly, one indeed finds that an exciton may behave effectively harmonic even in the presence 
of on-site anharmonicity, provided its delocalization length is sufficiently long. In fact, one can show 
by a perturbative expansion of the Hamiltonian Eq. (|56] |, treating the anharmonicity A as small, that the 
effective anharmonicity of an exciton scales like the so-called participation ratio lt20J : 
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Fig. 10. Pump-probe spectra of the C=0 mode of crystalline ACN as a function of temperature from 93 K to room 
temperature. Adapted from Ref. L20J . 



where 'Py are the expansion coefficients of the exciton in a site basis. The participation ratio is a com- 
monly used measure of the degree of delocalization, i.e. Lj^^ = 1 for ^ fully localized state and 
Y.j^'] ^ for a perfectly delocalized state. The exciton's anharmonicity is a direct measure of its 
degree of delocalization which, according to Fig.|9^, bottom, is large at 93 K. 

Interesting in this context is the temperature dependence of this effect (Fig.fTOli: At low temperatures, 
in the absence of disorder, the exciton is fully delocalized, and one can hardly bleach it. However, as 
temperature rises, disorder increases and the exciton starts to Anderson-localize. As a result, the exciton 
transition could then indeed be bleached because its effective anharmonicity increases (Fig. [TOl top). 
The onset of disorder-induced Anderson localization of the free exciton is observed at about the same 
temperature where the intensity of the self-trapped state disappears. This finding suggests that both 
localization mechanisms are connected. Disorder localization is caused by random variations of the 
diagonal elements of the excitonic coupling Hamiltonian Eq. ( |56] |. The amide I frequency is known 
to vary significantly with hydrogen bond distance (on the order of 20-30 cm^VA); hence, thermal 
excitation of phonons that give rise to random fluctuations of the hydrogen bond distances lead to 
disorder-induced localization of the free exciton. Since these are exactly the phonons which also mediate 
self-trapping at low temperatures, thermal excitation at the same time diminish the intensity of the 
self-localized state. Hence, it is the same nonlinear interaction, the variation of the amide I excitation 
energy with hydrogen bond distance, which gives rise to the two effects: (i) diagonal disorder and as 
a consequence, disorder-induced localization of the free exciton with rising temperature and (ii) self- 
localization at sufficiently low temperatures. The competition between Anderson (disorder) localization 
and nonlinear self-localization in the presence of disorder has been investigated in Ref. |42 1 : Those states 
that are self-localized for large enough nonlinear exciton-phonon coupling % are quite distinct from 
those that are Anderson (disorder) localized. Upon reducing the nonlinear coupling, the self-trapped 
states become less localized and eventually reach bifurcation points, below which they do not exist. 
Surprisingly, the states that are Anderson localized for vanishing nonlinear coupling become rapidly 
delocalized and unstable as nonlinear coupling is increased just slightly. 

Hamm and coworkers |20| concluded that the pump-probe experiments are fully consistent with 
earlier interpretations of conventional temperature dependent absorption spectroscopy ifTOI . namely that 
the higher frequency band of ACN relates to an free exciton whereas the lower frequency band relates to 
a self-localized state. Nevertheless, it turns out that the lifetime of either of these sates is short, 2 ps |20|, 
and as such, the original idea of Davydov l?) that vibrational self-trapping stabilizes the excitation, 
did probably not come true. A lifetime of 2 ps is a typical value for vibrational energy dissipation in 
condensed phase systems, and vibrational self-trapping is not capable to extend this lifetime by any 
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significant amount. 2 ps appears to be too short to trigger any serious conformational transition of a 
protein that could potentially be important in the biological function of the system. 

7.2 Comparison to a Numerically Exact Theory 

Although the discussion outlined above seemed to explain the experiment conclusively, we note again 
that the 'theory' of Ref. |20| artificially tried to 'cut out' the exciton part from the remainder. At the 
same time, it is clear that the free exciton is not an eigenstate of the full Hamiltonian Eq. ([BSl, which 
can be seen from simple insertion. Furthermore, also the experiments show that exciton and polaron 
are, of course, coupled to each other, leading to a spectral response of the latter when pumping the 
former (Fig.|9t, bottom). Unfortunately, none of the limiting regimes discussed in the previous sections 
apply to the parameter range that has been proposed for the ACN problem (% w 25 cm^ ' , co w50 cm^ ^ , 
y =4-10cm-i): 

- The DNLS (Sec. [3] and 15.2b is valid in the limit x/Ti(0 ^ °o, in which case the polaron binding 
energy Et, = —X^/Ti(x) would become large compared to the discreteness of the oscillation frequency 
h(0. However, for the given set of parameters for ACN we find x/h(0 = 0.5; by no means large! 
In fact. Fig. [8] shows that polaron solutions are expected in this regime for relatively large exciton 
couplings J, where the DNLS would strictly predict exciton solutions. 

- The small J limit seems to be better justified (Sec. 14.11 and 14.2b . Yet, this limit does not predict 
the existence of a free exciton, in disagreement with the common picture of the spectroscopy of 
ACN, and in disagreement with the interpretation of the experimental pump-probe results discussed 
above. According to the small / limit, one would rather expect a first excited state that lies exactly 
one phonon quantum above the ground state (Eq.l3Tb. However, a hypothetical phonon frequency of 
CO = 15 cm^^ (the spacing between two bands in Fig.|2]) would result in a completely different tem- 
perature dependence with a much lower turn-over between low and high temperature regime [,10|, 
and would furthermore contradict the experimental observation that it is indeed a co = 50 cm^^ 
phonon that couples strongly to the hydrogen bond 11912211 . 

The existence of the two solutions of the DNLS, an extended free exciton and an self-localized small 
polaron, was commonly used to rationalize the two peaks in the absorption spectrum of ACN (Fig.|2^). 
According to this interpretation, the lower frequency band relates to a self-trapped state with a temper- 
ature dependence given by the displaces oscillator picture 1 13 14 10| while the higher frequency band 
relates to the free exciton. Again, the free exciton is not an eigenstate of the full-quantum Hamiltonian 
Eq. ([T]-01i. Hence, we are in the bizarre situation that the DNLS apparently describes the experimental 
results correctly, knowing that we are in a parameter regime where it is not valid, while the formally 
more accurate full-quantum Hamiltonian does not. In light of this discussion, it might seem that the 
correct result of the DNLS is an unfortunate coincidence. 

In order to ultimately resolve this issue, Hamm et al. set out to calculate the absorption spectrum of 
ACN on a numerically exact level |36|, along the lines of Sec.|6l In doing so they were guided by the 
DNLS result that polaron and free exciton may coexist only in 3D, in which case barriers separate both 
solutions (see Sec. |6ll, and hence render an exciton solution meta-stable, even if it is slightly higher in 
energy. This is in accordance with the well known result from a simple quantum mechanical treatment of 
a square well potential in three dimensions, where a bound state appears only if the product Eqo^ exceeds 
a threshold value (with Eq the depth and a the width of the well). In contrast, an one-dimensional square 
well has at least one bound state for any value of Eqg^ |43 |. 

Hamm et al. were furthermore guided by the observation obtained from inspection of the 3D struc- 
ture of the ACN crystal (Fig. [T^) that exciton coupling along the hydrogen bond direction is negative 
(7|| > 0), while the coupling between hydrogen bond chains is expected to be positive (/x < 0). This is 
since exciton coupling is dominated by through-space transition-dipole coupling [44,45]; 

7 = -^^ + 3i^^i:^#L^, (58) 

r- ■ r- 

I] I J 

where /Ji and /jj are the C=0 transition dipoles of two sites (as vectors) and ry is a vector connecting 
them. The nearest neighbor distance r,j along the hydrogen-bonded chain is essentially the same as that 
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Fig. 11. Absorption spectrum of the Holstein in 3D with hai = 50 cm^', % = —25 cm^', y|| = 10 cm^' and 
= —10 cm^'. Stick spectra and spectra convoluted with a Lorentzian lineshape function with width 10 cm^' 
(FWHM) are shown. The eigenstates appear in clusters with properties of a polaron in its ground state (WPl), an 
exciton in its ground states (WP2), and a polaron in its first excited state (WP3, see Fig|12t. Adapted from Ref. [361 . 



between chains |46 |, which is why inter- and intra-chain couplings are expected to be of the same order, 
albeit with opposite signs. There is no reason to assume that inter-chain coupling can be neglected. 

Introducing anisotropy into the exciton coupling increases the parameter space considerably. In 
Ref. |[36| . an extensive search in parameter space was performed, staying within limits that seemed 
physically meaningful for the ACN problem, and looking for full-quantum solutions that predict co- 
existence of free exciton and polaron. Fig. [TT] shows the best result that could be obtained. The total 
oscillator strength is distributed over many eigenstates, nevertheless, the oscillator strength clusters in 
a way that two dominant bands remain when convoluting the stick spectrum with a Lorentzian line- 
shape function of width 10 cm^' (FWHM). The thus obtained spectrum resembles the experimental 
30 K spectrum reasonably well (Fig.|2l). The separation of the two bands (w 20 — 30 cm^') is somewhat 
larger than in the experiment, yet it clearly is smaller than one phonon quantum h(0. Furthermore, the 
higher frequency band shows some substructure in the simulation, similar to the experiment. 

We pause and note that we do not observe stationary eigenstates in a femtosecond pump-probe ex- 
periment (which are not resolved anyway in reality due to broadening mechanisms). Rather, a femtosec- 
ond laser pulse is spectrally broad (due to its shortness) and impulsively excites a coherent superposition 
state (wavepacket) of all eigenstates that are covered by the spectral width of the pump-pulse. The typ- 
ical width of the pulses used in Refs. I(20i2 1 J was 10 cm^', which is why it makes sense to discuss 
the properties of the wavepackets WPl, WP2 and WP3 (as indicated in Fig. [TTI). rather than that of 
individual eigenstates. 

A closer analysis shows that the eigenstates that cluster in the lower frequency band (WPl) resemble 
a ground state harmonic oscillator wavefunction with an origin displacement of about 1 (see Fig. fT2h). 
which is somewhat larger than the value predicted from the simple displaced oscillator picture with 
q — \/2x/(i) w 0.71. The eigenstates that cluster in the higher frequency band, on the other hand, can 
be classified in the following way: The lower energy portion of these states, clustering as a shoulder in 
the absorption spectrum (WP2), resembles ground state harmonic oscillator wavefunctions without any 
displacement, while the higher energy portion (WP3) resembles first excited state harmonic oscillator 
wavefunctions, again with origin shift. 

Fig- Eh shows the density of phonons around an exciton, calculated as 



d{j) = {^'\BlBob]bj\^') 



(59) 



Projected 
Wavefunction 



Phonon Density 




Fig. 12. (a) Projection of (t=0) wavepackets onto phonon coordinate q and (b) tlie corresponding plionon-densities 
d[j) = {'i'lBlBob'^jbjl^'). Wavepackets WPl, WP2 and WP3 are constructed from the embraced eigenstates as 
shown in Fig.[TT] The parameters in this simulation were: h(£> = 50 cm^', x = 25 cm^', = 10 cm^', = 
- 10 cm" 1 . Adapted from Ref. (36). 



which reveals essentially the same result: Wavepackets WPl and WP3 represent phonons dressed around 
an excitons, i.e. ground and first excited state of a polaron, whereas wavepacket WP2 has much smaller 
phonon contribution (albeit not zero). Hence, the results seem to be in reasonable agreement with ex- 
perimental observation. However, it should be noted that such free-exciton like solutions are found only 
in a very limited parameter regime. In particular, they are found only when exciton couplings with alter- 
nating signs of y|| and J± are introduced. The appearance of the double peak structure in the absorption 
spectrum of crystalline ACN thus seems to depend relatively critically on a balance of parameters, which 
might explain why vibrational self-trapping is observed so clearly only in ACN, while related crystals 
show, if at all, only minor effects ll22l . It is important to note that equivalent exciton-like solutions could 
not be found in ID, where the displaced oscillator picture (Sec. l4.2l l remains intact even for relatively 
large exciton couplings J. 

In Sec.|6]we have seen that exciton and polaron solutions may coexist in certain parameter regimes 
in the 3D case, as a result of a barrier that separates the two regions (Fig.|5]l. We find such a barrier only 
in the regime where the adiabatic DNLS becomes valid, i.e. for x/hoi ^ 1, while the transition between 
exciton and polaron is smooth in the quantum regime (Fig. [8]l. This can be rationalized as follows: In 
the quantum case (x/hoi < 1), the potential phonon displacement q is of the order of one or smaller, 
where the phonon displacement is defined in units of quantum mechanical delocalization length. In 
other words, quantum-mechanics would average out any barrier that occurs on length scales smaller 
than the delocalization length of the phonon wavefunction. In that sense, the coexistence of polaron and 
exciton solutions found in the numerically exact solution of the Holstein Hamiltonian (Fig.fTSl) is due to 
an effect that is different from the one discussed in Sec. |6l It is not a barrier on an effective (adiabatic) 
ground state surface (Fig. [T3b). but must be related to at least two adiabatic surfaces, as sketched in 
Fig. [T3b . A physical picture of these two surfaces, apart from the black-box result from the numerically 
exact diagonalization of the full Hamiltonian, still remains to be developed. Nevertheless, clearly it is the 
3D nature of the problem that is responsible for exciton-like solutions. In ID, the displaced oscillator 
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Fig. 13. Two mechanisms that might lead to coexistence of exciton and polaron solutions: (a) two regimes on one 
effective (adiabatic) potential energy surface separated by a barrier which, however, requires that the length scale 
of quantum mechanical delocalization of the wavefunction is smaller than that of the barrier, or (b) two effective 
(adiabatic) surfaces. 

picture (Sec. I4.2l i. which does not predict any exciton solutions, remains intact in a wide parameter 
range. 



8 Summary and Conclusions 

From its first appearance in the literature in 1983, the temperature dependent peak of the C=0 vibrations 
in acetanilide has been associated to the Davydov soliton idea. The latter assumes that when excitons 
are coupled strongly to acoustic phonons a soliton-like excitation forms that may propagate ballistically. 
Although Scott and collaborators used optical phonons for acetanilide, their selftrapping approach was 
similar to that of Davydov. The Davydov soliton idea has been very appealing from the very beginning 
since it asserts that nonlinear type modes play an important role in energy transfer in biology. Nonlinear- 
ity, of course, is not unknown to biology; on the contrary most functional motions of macromolecules 
are large amplitude ones and thus nonlinear. The theoretical as well as experimental fascination with 
acetanilide over a period of more than two decades results fundamentally from its similarity to peptide 
systems. If soliton-like objects exist in acetanilide then one may investigate their presence and function 
in proteins. 

On the theoretical side, the basic model of Scott for acetanilide involves the following elements: 
(i) The system is assumed to be one dimensional. This is because the C=0 dipoles along the hydrogen- 
bond axis seem to interact stronger among them than with the ones in the perpendicular directions, (ii) A 
bear vibrational C=0 excitation that tunnels slowly from bond to bond (/ = 5 cm^ ' ) couples to a single 
optical phonon (or to few phonons) of frequency (0 = 50 cm^ ' with coupling % = 25 cm^ ' . As a result a 
localized small polaron of the type of a DNLS discrete breather forms, (iii) The double absorption peak 
in the amide 1 region results from the polaron state at 1650 cm^' and the free exciton at 1665 cm^' (iv) 
The temperature dependence of the polaron state arises from the Frank-Condon factors connecting the 
unexcited phonon ground state with the displaced phonons in the excited polaron state. This theoretical 
Scott picture that is adiabatic in nature since it employs the Davydov methodology, was corroborated 
by Alexander and Krumhansl with the modification that acoustic phonons are also involved. 

In order to converge to the polaron picture, numerous alternative explanations had to be eliminated. 
The most serious objection to the polaron comes for the Fermi-resonance idea, viz. the accidental de- 
generacy of the amide 1 symmetric stretch mode with a combination band of an in-plane N-H bend 
deformation as well as a torsional mode 1,401 . Although Raman as well as neutron scattering experi- 
ments by Barthes and collaborators largely eliminated this possibility, the mere fact that upon complete 
deuteration the anomalous peak disappears, presents a yet unexplained feature that could speak in favor 
of a Fermi resonance 1 17 47 1. An alternative issue is related to the active presence of acoustic modes 
in selftrapping; although there is indication of their absence from incoherent neutron scattering, there is 
no direct measurement yet of the presumed optical modes that are excited in the processes of the C=0 
vibrations. While these two experimental issues seem to be still pending all other questions regarding 
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the role of low IR modes, presence of defects or possible hydrogen bond non-degeneracies have been 
conclusively addressed by the group of Barthes. 

While most of the absorption and scattering experiments were done in the eighties and nineties, 
the recent pump-probe experiments of Hamm and collaborators provide more direct information on 
the nature of the states involved. Through experiments in the acetanilide N-H band was possible to 
identify a 50 cm^' hydrogen bond related mode that assists selftrapping in this band. One then makes 
the natural assumption (compatible with the earlier work) that this phonon is responsible for the amide 
I selftrapping as well. Furthermore, in pump-probe experiments in the amide I region itself, they found 
that while the low frequency 1650 cm^' band corresponds indeed to a nonlinear mode, the 1665 cm^' 
one is certainly a linear one. As result, the former may be assigned to a selftrapped polaron state while 
the latter to a free exciton one, in a manner compatible with adiabatic theory of Scott. 

This adiabatic theoretical polaron scheme has been challenged recently by an exact diagonaliza- 
tion performed on the Holstein Hamiltonian as applied to the specific acetanilide problem. Using a 
numerically exact procedure Hamm et al. found that while the polaron state could be identified as the 
lowest energy state, the free exciton state did not appear in the calculated spectrum |36|. These numer- 
ics showed that there is no free exciton state neither in one nor in three dimensions for the "standard" 
acetanilide parameter values with hopping rate in the range 7 w 5 — 10 cm^'. In order to bypass this 
problem, Hamm et al. introduced different signs in the hopping rares depending on direction, viz. nega- 
tive overall rates along the hydrogen-bonded axis while positive along the plane perpendicular to it. The 
resulting spectrum is then compatible with the experiment within the standard parameter regime. 

The issue of the free exciton state has been discussed extensively in the present work. Using the 
semiclassical theory of Kalosakas et al., we showed that in three dimensions free exciton and polaron 
states are separated by a barrier while there is no such baiTier in one dimension |27|. Furthermore, for 
the standard acetanilide parameter regime the one dimensional semiclassical free exciton state has quite 
a short lifetime and is in general unstable. Thus, even semiclassically, the one-dimensional acetanilide 
picture needs to be revised. Using the double well adiabatic potential for free exciton and polaron and 
upon an ad-hoc requantization we find that the resulting two lowest states have polaron and free-exciton 
features respectively. Thus, within this intuitive, yet approximate, 3D picture we recapture in principle 
Scott's original explanation. The problem is that this result does not appear to agree quantitatively with 
the exact numerics mentioned above, unless different sign hopping terms are used in the latter. 

Cristaline acetanilide proved to be a good medium for experiment as well as theory for over twenty 
five years. During this period it permit the application of diverse experimental techniques as well as the 
testing of a large variety of theoretical ideas. It is clear that although many issues have been addressed 
and resolved, there are still many challenges ahead. In addition to the central dimensionality issue related 
to the free exciton, the absence of a temperature dependent peak from the fully deuterated acetanilide or 
the dispersion features of the modes that assist selftrapping need to be understood better. Furthermore, 
selftrapping in the N-H bond needs also to be addressed that is actually different compared to that in 
the C=0 bonds. Specifically, while a phonon progression that is absent in the later is seen in the former, 
it nevertheless terminates abruptly in a way dissimilar to what is expected from Eq. |3T] And while the 
C=0 overtones where treated by Scott through QDNLS, it would be interesting to analyze them directly 
through the Holstein model. It is clear that the exciting "acetanilide story" is far from over. 
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9 APPENDIX : Brief historical guide to the research in ACN 

Due to the numerous number of works published on the "Davydov soliton" as well as ACN, it is often 
hard to keep track of the various contributions and the significant developments in the "acetanilide 
story". We present here a partial list of contributions in a historical order that are specifically relevant to 
the acetanilide developments. 



Will be inserted by the editor 



27 



- 1973 Davydov introduces the idea that as a result to coupling with acoustic phonons a vibrational 
excitation may selftrap and become a non-topological soliton |3J. Careri finds experimentally an 
"anomalous peak" in acetanilide 16|. 

- 1982 Davydov publishes his book on "Biology and Quantum Mechanics" where he details his ap- 
proach to energy transfer in biomolecules 1231 . 

- 1983 Careri, Scott and collaborators present experimental and theoretical analysis of the C=0 band 
of crystalline acetanilide |7|. After eliminating possible alternative explanations they conclude that 
the temperature dependent ACN peak at 1650 cm^' is due to a "Davydov soliton" type mechanism. 

- 1984 Careri, Scott and collaborators present detailed experimental and theoretical details on the 
acetanilide "anomalous peak" |8|. Eilbeck, Lomdahl and Scott analyze selftrapping of the amide I 
mode through the Discrete Selftrapping Equation |9|. 

- 1985 Alexander treats theoretically the Davydov Hamiltonian using Holstein polaron techniques 1 1 3 1 . 
Blanchet and Fincher find experimentally that low IR as well as a band at 3250 cm^ ' have the same 
temperature dependence as that of the "anomalous peak" and suggest that a topological defect might 
be responsible for the anomalous amide structure |39|. Johnston and Swanson propose that that the 
anomalous peak is a result of Fermi resonance between the amide I mode and a combination mode 
between an N-H motion and a strongly temperature dependent low frequency phonon ll40l . Scott and 
collaborators introduce the quantum DST equation and explain higher frequency anomalous lines as 
overtone spectra of the C=0 band |41 1. 

- 1986 Takeno introduces a classical vibron model and studies the acetanilide temperature dependence 
in a way similar to that of localized modes in alkali halides |48 1. Alexander and Krumhansl provide 
an adiabatic treatment of the Davydov Hamiltonian and analyze the ACN peak using an acoustic 
phonon chain |14 |. 

- 1988 Barthes and coworkers study experimentally acetanilide and while confirm the presence of 
anharmonicity in the solid they report no evidence of acoustic phonons involved in selftrapping |24|. 

- 1989 Scott uses the theory of color centers and proposes a quantitative explanation for the tempera- 
ture dependence of the "anomalous" ACN peak ifTsll . 

- 1990 Austin and collaborators perform pump-probe experiments in acetanilide and propose that the 
temperature dependence of the anomalous peak is due to hydrogen bond non-degeneracy 1291 . 

- 1991 Barthes and collaborators show through inelastic neutron scattering experiments that ACN data 
are not consistent with a topological defect assumption 1 17|. They also claim that no indication of 
acoustic phonons are seen in the low IR spectra. Furthermore through Raman scattering experiments 
they report that while fully deuterated acetanilide shows no anomalous peak, a partially deuterated 
one does have one, downshifted however by aboutl50 cm^' from the main amide I peak. 

- 1992 Scott summarizes the "Davydov soliton" and the "acetanilide story" in a review article in 
Physics Reports |10|. 

- 1995 Barthes and collaborators use neutron scattering experiments and find no non-degeneracy in 
the hydrogen-bonded proton ll30l . 

- 1999 Kalosakas, Aubry and Tsironis investigate the semiclassical Holstein model and propose 
that phononic normal modes of the semiclassical polaron may be related to low frequency ACN 
modes |2749|. 

- 2000 Barthes and collaborators show that the low IR modes investigated earlier have no anomalous 
temperature dependence ll50l . 

- 2002 Hamm, Edler and Scott use pump-probe experiments and analyze the NH spectrum of ACN fT9l . 
They identify a phonon of approximately 50 cm^^ that seems to be involved in selftrapping of 
the exciton. Hamm and Edler publish pump-probe experiments for the C=0 band that identify the 
"anomalous" peak as a nonlinear state and the "normal" one as a linear free exciton state |20|. 

- 2006 Hamm and collaborator present a numerical diagonalization of the Holstein model in one and 
three dimensions and find inconsistencies between numerical results in one dimension and standard 
adiabatic anomalous peak explanation. They propose that the C=0 polaron are a result of the three 
dimensional nature of the ACN crystal as well as specific exciton dipole-dipole matrix elements l36l . 
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